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ABSTRACT 


Many  econometric  analyses  include  dependent  variables  which  are  constrained 
to  the  interval  between  zero  and  one.  Under  such  circumstances  simple 
regression  procedures  break  down.  The  simple  logit  model  is  extended  to 
the  multi-factor  case.  Two  alternative  models  are  defined  depending  on 
the  error  structure.  The  generalized  least  squares  approach  assumes  the 
share  specification  is  an  accurate  representation  of  the  underlying  input 
demand  structure.  The  multinomial  maximum  likelihood  model  treats  the 
dependent  variable  as  a  probability  with  a  multinomial  density.  Either 
model  provides  a  wide  array  of  functional  forms  while  maintaining  a 
straight-forward  estimating  procedure. 
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MULTINOMIAL  LOGIT:  A  LIMITED  DEPENDENT  VARIABLE  TECHNIQUE 

By 

MICHAEL  LEBLANC 
SEPTEMBER  1980 


I.  INTRODUCTION 

Many  econometric  analyses  include  dependent  variables  which  are  con¬ 
strained  to  the  interval  between  zero  and  one.  Under  these  conditions 
simple  regression  procedures  breakdown.  Predictions  may  fall  outside  the 
zero  one  interval  and,  because  of  the  grouped  nature  of  the  data,  error 
terms  are  heteroscedastic . 

The  objective  of  this  paper  is  to  outline  a  multinomial  logit  metho¬ 
dology  as  a  sensible  response  to  the  limited  dependent  variable  problem 
and  to  provide  a  more  robust  alternative  than  translog  cost  function 
approaches  to  estimating  derived  demand. 

The  first  section  of  the  paper  describes  the  limited  dependent  variable 
problem  and  the  two  approaches,  probit  and  logit,  typically  used  to  resolve 
this  difficulty.  The  simple  logit  model  is  extended  to  the  multi-factor 
case  and  discussed  in  a  probability  model  context.  Two  alternative  models 
are  suggested  and  estimating  methodology  for  each  is  discussed.  Finally, 
the  implications  for  normalizing  the  models  to  facilitate  estimation  is 
addressed. 

II.  LIMITED  DEPENDENT  VARIABLE  MODELS 

Models  with  limited  dependent  variables  have  been  used  in  many  fields 
including  biology,  psychology,  and  economics.  The  earliest  work  was  Berkson' 
(1944,  1951,  1953,  1955)  applications  to  bioassay  problems.  The  initial 
analysis  in  economics  was  Tobins  (1958)  automobile  demand  study.  More 
recent  applications  using  multinomial  logit  methodology  are  Tyrrell's 
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(1979)  analysis  of  household  expenditures,  Baughman  and  Joskow’ s  (1975) 
study  of  appliance  choice,  and  Joskow  and  Mishkin’s  (1977)  electric 
utility  fuel  choice  model.  In  the  case  of  a  two  variable  probability 
choice  model  the  simplest  expression  is 

(1)  yt  =  a  +  Bxt  +  et 

where  Yt  is  a  proxy  for  unobserved  the  probability  Pt,  xt  is  an 
explanatory  variable,  and  et  is  the  error  term.  The  analyst’s  interest 
is  the  conditional  probability  of  decision  yt,  purchasing  capital 
equipment,  for  example,  given  explanatory  value  xt,  the  equipment  price. 

It  is  not  possible  to  know  whether  a  particular  firm  will  purchase  equipment 
when  the  equipment’s  price  is  xt.  The  objective  of  the  analysis  is  to 
determine  the  probability  of  purchase.  This  implies  a  dependent  variable 
confined  to  the  interval  from  zero  to  one. 

Simple  linear  regression  equation  (1)  may  result  in  prediction 
problems  (See  Figure  1). 


Figure  1.  SIMPLE  REGRESSION 

While  yt  is  constrained  between  zero  and  one  a  +  8xt  is  unconstrained. 
It  is  quite  possible,  especially  in  a  forecast  mode,  for  predicted  y 
to  fall  outside  the  zero-one  bounds.  In  addition,  since  the  residual 
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variance  depends  on  the  unknown  probability  value  Pt  the  variances  are 
heteroscedastic. 


To  mitigate  the  prediction  problems  a  monotonic  transformation  is 
applied  to  expand  the  probability  range  from  negative  to  positive 
infinity  (See  Figure  2). 
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FIGURE  2.  MONOTONIC  TRANFORMATION 
OF  A  PROBATILITY  TO  THE 
RANGE  (-  00  ,°°  ). 


The  probability  P  is  measured  on  the  vertical  axis  and  its  transformation 
along  the  horizontal.  As  P  increases  from  zero  to  one,  its  transformation 
increases  from  negative  to  positive  infinity.  The  prediction  problem  of 
using  a  limited  dependent  variable  is  therefore  avoided.  There  are  many 
different  transformations  that  may  be  applied.  The  most  common  transfor¬ 
mations  utilize  either  the  normal  distribution  or  logistic  function. 

Probit 

By  making  a  specific  assumption  about  the  distribution  of  the  error 
term  the  parameters  of  equation  (1)  can  be  estimated.  If  the  culmulative 
distribution  function  of  a  standardized  normal  variate  is  used  the 
technique  is  termed  probit  analysis.  The  critical  level  of  x  triggering 
a  purchase  decision  varies  from  observation  to  observation.  By  assuming 
the  distribution  of  these  critical  levels  is  normal,  the  conditional 
probability,  P(x),  can  be  expressed.  The  dependent  variable  y  is 
implicitly  defined  by 

a  +  6  x 

1/2 

(2)  F(y )  =  P(x)  =  ( 1/ 2tt  )  /  exp  ((-e/2)/2)de 
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where  P  is  the  probability.  The  transformation  used  is 
-1 

(3)  y  =  F  (P)  =  ct  +  3x. 

While  circumventing  the  prediction  difficulty  the  transformation  has  not 
mitigated  the  heteroscedasticity  problem. 

The  transformed  model's  parameters  can  be  estimated  by  maximum  likeli¬ 
hood  methods.  Let  x  xn  be  observations  for  the  independent  variable. 

By  arranging  the  observations  where  the  first  n*  firms  are  purchasers 
and  the  last  n  -  n*  are  not,  the  logarithmic  likelihood  function  is 
n*  n 

(4)  InL  =  ^  In  P(x^)  +  2  ln(l  -  P(xt)) 

t=l  t=n*+l 

where 

1/2  2  2 

(5)  (1/2 Trcr )  exp(-(yt  -  a  -  6xt)  /2 o  )) 

The  maximum  likelihood  estimates  are  obtained  by  maximizing  InL  with 
respect  to  a,  3,  and  a2.  The  solution  may  be  obtained  using  non-linear 
techniques. 

Logit 

Another  transformation  utilizes  the  logistic  function.  In  a  two  state 
decision  framework  the  logistic  is  written 

(6)  Yt  =  1/(1  -  exp  ( g ( x )  +  et)) 

where  Yt  is  defined  above  and  g(xt)  is  an  undefined  function.  By  writing 
the  odds  in  favor  of  equipment  purchase  as  the  ratio, 


(7)  P/(l-P)  =  exp(-  (g(x)  +  e)) 
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where  P  is  the  purchase  probability,  the  range  of  the  dependent  variable  is 
expanded  from  negative  to  positive  infinity.  The  estimated  form  is  obtained 
by  taking  the  logarithm  of  equation  (7) 

(8)  In  [P/(l-P)]  =  g(x)  +  e. 

The  logarithm  of  the  probability  ratio  is  a  monotonically  increasing  function 
of  the  probability  P  varing  between  negative  and  positive  infinity.  It  is 
numerically  equal  to  the  logit  of  the  complementary  event  but  opposite  in 
sign 

(9)  In  (P/( 1-P) )  =  -  In  ((1-P)/P). 

To  estimate  the  logit  function  the  error  term  is  assumed  to  be  generated  by 
the  binomial  distribution.  The  odds  or  chance  of  success,  equipment 
purchase,  for  any  observation  is  analogous  to  sampling  with  replacement. 

The  binomial  distribution  is 

(n\  w  n-w 
w /  P  (1-P) 

where  n  is  the  total  number  of  dollars,  P  is  the  probability  of  success, 
and  w  is  the  number  of  samples  taken,  share  of  equipment  purchasers  in 
the  sample.  If  by  specifying  P^  as  a  logistic  the  likelihood  function 
is 

n  rn\  w 

(11)  L  =  n  lw/ (1/(1  +  exp  (-g(xt))))  . 

t=l  n-w 

(1  -(exp  g(xt)/(l  +  exp  (-g(xt))))) 

The  density  function  f(w)  is  interpreted  as  the  probability  of  choosing 
w  positive  responses  from  a  sample  of  n. 
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The  ability  of  the  logistic  function  to  constrain  probability  esti¬ 
mates  to  the  range  between  zero  and  one  is  an  important  reason  for  the 
popularity  of  this  model.  In  addition,  the  logit  model  is  easy  to  estimate, 
especially  when  the  logarithmic  transformation  is  linear  in  the  parameters. 
Furthermore,  the  underlying  probability  function  gives  a  close  approxima¬ 
tion  to  the  cumulative  normal  probability  function. 

III.  MULTINOMIAL  LOGIT 

While  the  two  state  logit  decision  model  provides  an  appropriate 
beginning,  its  usefulness  is  limited.  Henri  Theil  (1969)  first  proposed 
a  multinomial  extension  of  the  linear  logit  model  to  analyze  consumer 
budget  allocations.  The  multinomial  logit  is  a  flexible  approach  and 
for  many  applications  provides  a  useful  marriage  between  theoritically 
founded  models  and  more  ad  hoc  prediction  oriented  models.  Unlike  other 
share  approaches  for  production  or  consumer  behavior  the  multinomial 
logit  does  not  require  expenditures  to  be  optimal  or  the  result  of  a 
decision  process  having  a  unique  mathematical  representation.  Flexability 
is  manifested  in  the  ability  to  incorporate  many  factors  influencing 
expenditures.  Since  a  major  objective  of  this  paper  is  to  articulate  an 
alternative  to  translog  cost  functions  the  following  discussion  is  couched 
in  terms  of  derived  input  demand. 

The  multinomial  logit  model  indicates  producer  expenditures  on  factor 
inputs.  A  set  of  N  independent  conditional  probabilities  are  assumed  to 
jointly  determine  the  allocation  of  expenditures  into  N  categories.  In 
addition,  the  probabilities  are  assumed  to  have  the  logistic  structure  and 
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are  conditional  on  input  prices  and  other  independent  variables.  The 
probabilities  are  written 

n 

(12)  v±  =  expCg^CP ,x) )  /  i  exp(gk(P,x)) 

k=l 

where  n  is  the  conditional  probability  of  a  dollar  being  allocated  to  input 
i,  P's  are  input  prices,  x's  are  other  independent  variables,  and  g-[(P,x)'s 
are  generalized  functions.  No  specific  restrictions  on  g^(P,x)  are 
necessary  to  apply  the  multinomial  logit.  However,  to  aid  exposition 
the  functions  g£(P,x)  are  assumed  linear  in  their  parameters  and  invariant  in 
structure  between  equations  except  for  inter-equation  parameter  variability. 
The  functions  g^(P,x)  are  specified  as 

(13)  gi  =  6iih^(Pi)  +  Bi2h2(p2)  +  •••  +  + 

Yilgl(xi)  +  Yi282(x2)  +  •••  +  YimSmCxm) 

(i  =  1,2,...,  N). 


Since  is  unobservable  the  model  is  made  operational  by  using  input 
expenditures  or  cost  shares  as  a  proxy.  Equation  (12)  becomes 

N 

(14)  SjL  =  PjQ-j/M  =  exp(gi(P,x))/  E  exp(gi(P,x))  (i  =  1,2,...,  N) 

k=l 


where 


N 

(15)  M  =  £  PiQi. 

i=l 


To  derive  the  sensitivity  of  the  ith  input  share  to  changes  in  each 
independent  variable  the  slope  of  the  ith  budget  share  with  respect  to  any 


variable  can  be  determined  as 
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N 

(16)  S±/dy  =  SiOgi/ay  -  I  Si  ZgJdy)  (i  =  1,2 . N), 

j-l 

Caution  is  warranted  when  using  estimated  parameters  from  the  multinomial 
logit  model  to  derive  demand  price  sensitivity.  The  logit  approach  is  used 
to  estimate  parameters  characterizing  the  underlying  probability  function. 
Therefore,  using  the  parameters  to  characterize  price  sensitivity  presumes 
input  demand  allocations  are  distributed  in  the  same  manner  as  the 
corresponding  probabilities. 

Elasticity 

Two  alternative  derivations  for  price  elasticities  are  possible. 

First,  total  expenditures  are  assumed  invariant  to  input  price  changes. 

This  is  the  production  analogue  to  the  budget  allocation  problem  of 
consumer  theory.  The  logit  model  is  converted  to  a  system  of  input  demand 
equations 

N 

(17)  Qi  =  MSi/Pi  =  exp  (f±  +  InM  -  InPi)/  Z  exp(fj)  (i  =  1,2,...,  N). 

j-l 

Given  the  "budget"  assumption,  the  own-price  elasticities,  Ej^,  cross-price 
elasticities,  E^j,  and  elasticities  with  respect  to  all  non-price  independent 
variables  are  derived  for  the  ith  good 

N 

(18)  Ei±  =  Pi(3fi/3Pi  -  z  Si  af-t/aPi)  -  1  (i  -  1,2,...,  N) 

j-l 

N 

(19)  Eij  =  PjOfi/BPj  -  Z  Sj3fi/3Pj)  (i  =  j  =  1,2,...,  N) 

j-l 

N 

(20)  Elx  «  xk(3fi/3xk  -  I  S j 3f j/9xk)  (k  =  1,2,...,  M). 

k  j-l 
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The  own-price  elasticity  approaches  negative  one  and  the  cross¬ 
price  elasticity  approaches  zero  as  the  ith  input  share  approaches  one  or 
the  derivatives  3 f /  3P^  are  invariant  across  shares.  While  the  elasti¬ 
cities  approach  zero  or  one  there  are  no  a  priori  restrictions  on 
either  the  magnitude  or  the  sign  of  the  elasticities. 

If  we  assume,  however,  the  total  expenditure,  M,  is  not  invariant  to 
changes  in  input  prices  and  therefore  3M/3Pj  ^  o  then  the  own-price  and 
cross-price  elasticities  are 

N 

(21)  Eij*  -  wA  +  PiOfi/gPi  -  Z  SjSfj/SPi)  -  1  (i  =  1,2,...,  N) 

j=l 

N 

(22)  Ei:j*  =  Wi  +  Pj(3f1/3Pj  -  Z  SjSfj/SPj)  (i  =  1 ,2 , . . .  ,  N) . 

j=l 

In  either  case  the  calculation  of  Eix^  remains  the  same.  If  the  own- 
price  elasticities  are  negative,  the  effect  of  input  price  increases  on 
costs  is  to  reduce  price  sensitivity  by  reducing  total  output. 

IV.  ESTIMATING  THE  MULTINOMIAL  LOGIT 

Specifying  the  multinomial  logit  input  expenditure  model  requires 
selecting  a  functional  form  for  f^.  The  function  f^  includes  price  and 
output  effects  and  all  other  production  or  producer  characteristics  con¬ 
sidered  relevant.  Assuming  f ^  is  linear  in  parameters  facilitates 
linear  estimation  techniques. 

Two  stochastic  specifications  of  the  model  are  possible.  First,  the 
share  equation 

N 

(23)  Si  =  PiQi/M  =  exp(fi)/  Z  exp(fk)  (i  =  1,2,...,  N) 

k=l 
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can  be  directly  interpreted  as  a  accurate  characterization  of  producer  expen- 
ditures  for  inputs.  The  stochastic  model  structure  is  generated  by  appending 
error  terms  to  each  function  fj..  The  input  share  equation.  ,  is 
transformed  to  an  estimable  form  by  normalizing  with  share  N  and  taking 
logarithms.  When  f^  is  linear  in  parameters  estimated  equations  are 

(24)  ln(Si/SN)  =  (  %n>  +  (B±j- 

+  (BNi  “  6NN>h(pN)  +  (Yoi  ■  yon) 

+  (Yli  -  YiN)g(xi)  +  ...  +  (YNi  -  Ywi)g  (xH) 

+  (£i  -  eN)  (i  =  1,2,... ,  N). 

Clearly,  the  presence  of  in  each  equation  generates  non-zero 
covariance  between  equation  error  terms.  Estimating  each  equation 
separately  by  ordinary  least  squares  leads  to  inefficient  parameter 
estimates.  If  each  error,  ,  is  assumed  homoscedastic ,  a  generalized 
least  squares  estimator  can  be  applied  to  improve  estimator  efficiency 
relative  to  an  ordinary  least  squares  procedure. 

Zellner’s  (1962)  genralized  least  squares  technique  for  "Seemingly 
Unrelated  Regressions"  can  be  used  to  improve  the  efficiency  of  the  esti¬ 
mators.  If  the  same  set  of  regressors  is  utilized  in  all  N-l  equations 
then  Zellner's  procedure  and  ordinary  least  squares  gives  identical  coeffi¬ 
cient  estimates.  Estimates  of  the  standard  errors,  however,  differ. 

Parameters  estimated  by  the  Zellner  generalized  least  square  procedure 
are  not  invariant  to  the  choice  of  common  denominator  budget  share  when 
an  estimated  variance-covariance  matrix  is  employed.  Barten  (1969), 
has  shown  maximum  liklihood  estimates  are  invariant  to  which  equation 
is  deleted  and  Kementa  and  Gilbert  (1968)  and  Dhrymes  (1970)  show 
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iterating  the  Zellner  procedure  leads  to  maximum  likelihood  estimates. 
Therefore,  iterating  the  Zellner  procedure  is  a  computationally  efficient 
means  of  deriving  maximum  likelihood  estimates.  In  general,  the  properties 
of  maximum  likelihood  estimators  only  hold  asymptotically.  However,  Kmenta 
and  Gilbert  (1968)  suggest  that  most  of  the  maximum  likelihood  estimator's 
asymptotic  properties  tend  to  be  present  in  small  samples. 

The  second  interpretation  of  the  multinomial  logit  is  to  assume  S.^  is 
a  proxy  variable  for  the  unobservable  probability  ^  generated  by  the 
multinomial  distribution 

X1  x2  XN 

(25)  f(xi,x2,  ...,  xn)  /  (M!/(xi!x2!  ...  x^!))  Pi  P2  ...  ?N 
where 

N 

(26)  Z  Pi  =  1. 
i=l 

The  multinomial  distribution  characterizes  the  results  of  sampling  with 
replacement  from  a  population  classified  into  N  >  2  categories.*  The  sampling 
procedure  can  be  viewed  in  several  ways  corresponding  to  different  units 
of  measure.  The  implication  of  different  measurement  unit  is  to  alter  the 
weights  in  the  likelihood  function.  However,  the  maximum  of  each  likelihood 
function  is  satisfied  by  the  same  set  of  probability  estimates. 

The  multinomial  expression  utilized  for  the  input  demand  expenditure 
system  allocates  each  percentage  of  total  expenditures  independently.  If 
shares  are  multiplied  by  100  to  facilitate  the  use  of  factorial  notation 
then  the  multinomial  distribution  is 


*See  Fellner ,  W. ,  "An  Introduction  to  Probability  Theory  and  Its  Appli¬ 
cation,"  (John  Wiley  and  Sons,  1968)  p.  167. 
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lOOSi  100SN 

(27)  f ( x)  =  (100!/(100S1!100S2!  ...100Sn!))  77 1  ...*N 

for  period  t. 

Replacing  77  .j_  with  the  appropriate  logit  expression  gives 

(28)  f ( x )  =  (100!/(100S1!100S2!  ...100Sn!)). 

N  lOOSi  N  100SN 

(exp(f1)/  ^  exp(f1))  ...  (exp(fN)/Z  exp  (fN)) 
i=l  i=l 


The  joint  multinomial  distribution  is  used  to  generate  a  likelihood 
function  of  the  form 


(29) 


L  =  At  Z  Z  (exp(f.)/  £  exp  (f,)) 
t=i  i-i  k-i 


N 


N 


100S 


It 


where 


(30)  AT  =  Z  (100 ! /100Sit !  ...  100SNt!). 

t=l 

Since  the  maximum  likelihood  estimator  is  invariant  to  monotonic 
transformations,  the  likelihood  function  is  simplified  by  writing  it  in  its 
logarithmic  form 

T  N  N 

(31)  In  L  =  Constant  +  100  Z  (  Z  S^tfj_t  -  ln(  £  exp(f^t))) 

t=l  i=l  i=l 

The  log-likelihood  function  is  maximized  where  the  first  order  derivatives 
with  respect  to  the  parameters  and  the  matrix  of  second  order  derivatives 
is  negative  definite.  The  first  order  conditons  are 


BlnL/SBik  =  £  xit(J>it  -  ^^t)  =  0 

t=l 


(i  =  1,2,... ,  N) 
(k  =  1,2 . K). 


(32) 
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The  second  derivatives  are  therefore 

(33)  3  WBlk  '  xktxkt'nit(nit  -  1) 

t=l 

where  both  parameters,  8ik  and  8^^ appear  in  the  same  equation  and 

T 

(34)  9 2 InL/  98ik  36 ik*  =  1  xktxkt '  Hitn-jt 

t=l 

and  each  parameter  is  from  a  different  equation.  Since  the  first-order 
conditions  are  non-linear  an  iterative  search  procedure,  Newton-Raphson 
or  Davidon-Fletcher-Powell  for  example,  must  be  employed  to  solve  for 
the  maximum  likelihood  parameter  estimates. 

There  is  an  indeterminancy  in  the  equation  system  arising  from  the 
constraint 


N 

(35)  Z  S±  »  1. 

i=l 


The  share  equations  are  invariant  to  the  addition  of  the  same 
expression,  lnz  for  example,  to  every  fi  since 

N 

(36)  exp(g4  +  6lnz)/  Z  exp(f^  +  81nz)  = 

3=1 

N 

exp(gi)exp(6lnz)/(  Z  exp(gj)  (exp(6lnz)) 

3=1 

N 

=  exp(gi)/  Z  exp(gj)  =  Si  (i  -  1>2,  ...  N). 

3=1 

This  indeterminancy  causes  the  matrix  of  second  order  partial  deriva¬ 
tives  in  the  multinomial  maximum  likelihood  procedure  to  equal  zero. 

The  singularity  problem  is  avoided  by  normalizing  the  N  parameters  for 
a  particular  variable  in  the  N  functions  (Tyrrell,  1979).  The  normali— 
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zation  imposed  does  not  affect  the  computed  elasticities  or  the  predicted 
shares.  The  most  straight  forward  approach  is  to  divide  the  N-l  equations 
by  the  Nth  equation.  The  equations  are 

(37)  Sj/SN  =  exp( f ±  -  fN)  (i  =  1,2,  ...,  N-l). 

In  the  estimated  logarthmic  form  the  equations  are 

(38)  In  ( SjL / SN )  =  (B0i  “  B0n)  +  (B  li  "  iN^l^l^  + 

(6Ni  *  6NN)hN(pN)  +  (Y  oi  "Y0n)  +  (Y  li  “Yin) 
gl(xX)  +  ...  +  (YNi  ”Ynh)8m^xm)  (i  =  •••»  N). 

The  indeterminancy  caused  by  the  adding-up  constraint  prevents  the 
calculation  of  the  individual  B^k's  but  allows  the  derivation  of  both  the 
predicted  probabilities  and  the  elasticities.  To  facilitate  discussion 
assume  only  price  variables  are  included  in  the  estimated  model.  Therefore, 
hf(Pi)  equals  Pj_  where  P^  is  the  ith  price.  If  the  share  system  is 
normalized  using  the  Nth  share,  then  the  estimated  shares  are 

A  A  A 

(39)  In  S^*  =  Bio  +  B-qPi  +  ...  +  B^n  Pjg  (i=l,2,  ...,  N-l) 

where  S^*  is  equal  to  the  predicted  ratio  of  the  ith  to  the  Nth  share  and 

A 

Bij  is  equal  to  the  estimated  coefficient  for  the  ith  share  and  jth  variable. 
The  predicted  Nth  allocation  is  derived  by 

N-l 

(40)  SN  =  1/(1  +  E  exp(lnSi*)Sfl)) 

j=l 


since 

N-l 

(41)  1  =  1/(SN  +  (  £  exp(lnSi*)SN). 

j-1 


(«< 

. 
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The  other  N-l  allocations  are  easily  calculated  by 
(42)  =  SN  exp(lnS^*). 


To  derive  price  elasticities  the  Nth  elasticity  must  necessarily  be 
derived  first.  The  predicted  elasticity  for  the  Nth  allocation  with 
respect  to  the  ith  price  is 


*  N-l  A 

(43)  nNi  «  -  P±  E  Sj  df j*/3Pi  +  ( -1  if  i=N 

j=l  1  0  otherwise 

( i=l ,  2 , . . . ,  N). 

The  ability  to  estimate  the  Nth  rests  on  the  equivalence  of 


(44) 


N-l 

-  P.(  l  Si9fj*/9Pi)  =  PiOfj/9Pi 


N 

-  z  s-af./sPj^) 

i=l  J 


where  the  right-hand  side  of  equation  (44)  is  derived  from  the  unnormalized 
price  elasticity  formulation  given  in  equations  (18)  (19),  and  (20).  The 
two  components  in  equation  (44)  are  equivalent  since 
N-l  N-l 

(45)  -  P4(  Z  S1  afWaPi)  =  -  PA(  Z  Si  (8  i±  -  6Ni) ) 

J-l  j-1 


N-l 

=  -  pi(  2  sj  $ji  ~  ^Ni 

J-l 


N-l 
2  Si) 
J-l 


N-l 

“  -  pi(  Z  Sji  “  eNi(1"SN)) 

J-l 

N-l 

=  -  Pi(  E  Sj  Bji  -  8ni  +  Bni  Sn) 

J-l 

N 

“  pi^Ni  "  z  eji  sj)* 

J-l 


hji  =  rifli  +  Pj:9f  j*/3Pi 


(46) 


( i— 1 , 2 ,  . ..,  N— 1). 
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Therefore,  although  values  cannot  be  estimated  for  each  Bji  coefficient 
predicted  shares  and  elasticities  can  be  derived.  The  only  restriction  is 
that  in  each  case,  shares  and  elasticities,  the  Nth  value  is  calculated 
first. 

V.  CONCLUSION 

Limited  dependent  variable  estimating  techniques  are  receiving 
increased  attention  for  consumer  budget  and  derived  demand  input  studies. 
The  multinomial  logit  model  provides  a  flexible  alternative  to  the 
more  popular  translog  expenditure  approach  by  allowing  a  wide  range  of 
explanatory  variables  while  constraining  predictions  to  the  interval 
between  zero  and  one.  Two  general  models  are  defined  depending  on  the 
error  structure.  The  generalized  least  squares  approach  assumes  the 
share  specification  is  an  accurate  representation  of  the  underlying 
input  demand  structure.  The  multinomial  maximum  likelihood  model 
treats  the  dependent  variable  as  a  probability  with  a  multinomial 
density.  Either  model  provides  a  wide  array  of  functional  forms  while 
maintaining  a  straight-forward  estimating  procedure. 


I 
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